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SUMMARY 

A computer  program  has  been  written  to  determine  the  local  flow  field 
about  an  airplane  due  to  wing-body  effects.  The  results  are  in  terms 
of  local  angles  of  attack  and  sideslip  at  specified  points.  A simple 
geometric  model  is  used  with  body  effects  determined  from  slender  body 
theory,  circulation  effects  (including  wing-body  interference)  calculated 
using  a series  of  line  vortices,  and  wing  thickness  effects  calculated 
using  thin  wing  theory.  The  various  geometric  parameters  plus  the  flight 
condition  (in  terms  of  Mach  number,  lift  coefficient,  and  angle  of  attack) 
are  input.  (Span  and  chord  loading  may  also  be  varied  if  desired.) 
Compressibility  effects  are  accounted  for  but  are  limited  to  the  subsonic 
(below  critical)  Mach  number  range. 

A description  of  the  program,  including  the  assumptions  made  and 
model  used,  and  of  the  use  of  the  program  is  contained  in  this  report. 
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INTRODUCTION 

For  store  separation  studies  or  other  purposes,  it  is  often  desirable 
to  know  the  local  flow  field  about  an  airplane  (in  terms  of  flow  angularity). 
Wind  tunnel  tests,  although  best,  are  expensive  and  complicated  and  thus 
not  practical  in  many  cases.  The  alternative  is  some  form  of  numerical 
calculation.  Numerical  computations,  in  attempts  to  become  "exact",  have 
a history  of  becoming  very  complicated  without  reproducing  the  actual  flow, 
because  assumptions  and  modeling  still  have  to  be  made.  In  cognizance  of 
all  this,  a computer  program  was  developed  using  a relatively  simple  model 
for  determining  the  local. flow  field  about  an  airplane  due  to  wing  and 
body  effects.  The  program,  and  the  way  to  use  it,  are  described  in  this 
report. 

s,  •;  . *'  -t  > 'V  i t * 

PROGRAM  DESCRIPTION 

A computer  program  (PROGRAM  WBFLFD)  has  been  written  (for  a CDC  6600 
Computer)  to  obtain  rough  approximation  of  the  local  airplane  flow  field  due  to 
wing-body  effects.  The  results  are  in  terms  of  local  angles  of  attack  and 
sideslip  at  given  (inputted)  points.  Compressibility  effects  are  included 
but  are  limited  to  the  subsonic  (below  critical)  Mach  number  range. 

The  wing-body  geometric  model  used  is  a straight-tapered  wing  with 
no  dihedral  and  an  axisymmetric  body  which  may  be  at  any  height  or  incidence 
relative  to  one  another. 

The  model  used  to  generate  the  riow  field  is  as  follows:  The  local 
flow  perturbations  are  considered  to  arise  from  three  sources:  circulation 
effects  (mostly  from  the  *•*•}. but  inclvli  ig  effects  of  wing-body  interference), 
body  effects,  and  wing  thiw„as*  effects  The  circulation  effects  are 
calculated  by  summing  the  fie*  -lor  it  i ef  induced  by  a series  of  line  vortices. 
The  relative  strengths  of  -;'•*«!  * vortices  are  determined  such  that  desired 
spanwlse  and  chordwise  c trib^Mona  are  approximately  reproduced 

(including  the  effects  oi*  /ing  body  Interference).  The  absolute  strengths 
of  the  vortices  are  det iras  . 7ach  that  a specified  lift  coefficient  is 
produced.  Figure  1 shows  the  arrangement  of  the  line  vortices  on  the  wing 
and  trailing  behind  the  wing.  Note  that  the  theoretical  wing  (extending  to 
body  centerline)  is  used  and  that  there  are  ten  spanwise  (in  each  semispan) 
divisions,  between  which  vortex  strength  is  varied  and  trailing  vortices  y 
leave  the  wing,  and  three  chordwise  divisions.  The  spanwise  lift  distribution 
(in  terms  of  circulation  strength  T at  a given  spanwise  location  y)  is 
determined  as  follows:  For  an  elliptical  lift  distribution, 

r - r0  /l-(y/s)2 

where 


4 

r„  - - c 


( ) 
uv  ' 2 ; 


In  these  equations  s is  the  semispan,  Cl  is  the  lift  coefficient,  Cflv  is 

the  nvcrngu  chord,  and  U is  the  flight  velocity.  (For  the  present  application, 
U Is  equat  to  1.0  so  the  vertical  and  sideways  perturbation  velocity 
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components  which  are  calculated  are  equal  to  the  upwash  and  sldewash  angles.) 
Wing-body  Interference  effects  cause  a change  in  the  spanwlse  lift  distribution 
This  Is  modeled  in  the  program  as  shown  In  Figure  2.  The  constants  a 
and  F shown  there  are  determined  such  that  the  net  area  under  the  curve 
stays  constant  (that  Is,  the  same  as  In  the  elliptical  distribution),  which 
means  the  total  lift  stays  the  same,  and  also  such  that  the  ratio  of  exposed 
wing  lift  In  the  presence  of  the  body  to  body  lift  due  to  the  wing  is  In 
agreement  with  an  accepted  procedure  for  calculating  this.  (See,  for  example, 
reference  (a).)  The  spanwlse  lift  distribution  can  also  vary  due  to  such 
factors  as  twist,  sweep,  aspect  ratio,  taper  ratio,  and  Mach  number,  which 
could  cause  either  an  Inboard  or  an  outboard  shift  in  the  distribution 
(without  changing  the  total  lift).  Although  usually  the  effect  Is  small, 
this  can  be  accounted  for  In  the  program  by  inputting  a factor  E such  that 
at  any  spanwlse  location 


3.82(y/s)  ■+  0.955,  if  y/s  <0.5 
sin  2ir(y/s)  , if  y/s  >0.5 


The  variable  c.  Is  the  local  lift  coefficient  and  c Is  the  local  chord. 

• 1 
The  values  of  the  circulation  at  ten  equally  spaced  spanwlse  locations  (In 

each  semispan)  are  used  In  the  program  calculations.  (Some  charts  to 
estimate  span  loading  changes  are  contained  in  reference  (b) . ) 


In  the  chordwise  direction,  the  program  nominally  assumes  a flat  plate 
loading  distribution,  with  line  vortices  located  respectively  at  0.135c, 
0.25c,  and  0.63c  having  relative  magnitudes  of  0.503,  0.345,  and  0.152 
(sum  equals  1.0).  The  exact  distribution  Is  usually  not  critical,  but  may 
be  varied  In  the  program  by  adding  Increments  to  the  first  two  relative 
magnitude  values  given  above  (with  the  third  value  automatically  adjusted 
so  that  the  sum  remains  1.0). 


. Flow  perturbations  due  to  body  effects  are  calculated  In  the  program 
using  slender  body  theory.  The  equations  are: 


w and  v are  the  perturbation  velocity  components  In  the  vertical  (z)  and 
sideways  (y)  directions  respectively,  R(x)  is  the  body  radius  at  the 

x-locatlon  of  the  Inputted  point  In  question,  R^(x)  is  the  slope  of  the  bod; 
at  the  x-locatlon  of  the  Inputted  point  in  question,  r is  the  distance  from 
the  point  In  question  to  the  body  centerline, 
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where  z_  is  the  z location  of  the  fuselage  (or  body)  centerline  and  y is 
measured  from  the  airplane  plane  of  symmetry,  and  a_  is  the  effective  angle 

of  attack  of  the  body. 

• #\'. • • 5 • < ,u  : ; ,;i  - , --t  . 

Flow  perturbations  due  to  wing  thickness  are  generally  small  relative 
to  the  other  contributions.  Thus,  the  exact  thlcknec*  distribution  is  not 
critical  and  an  assumed  thickness  distribution  (having  a maximum  thickness 
located  at  0.368c)  has  been  used  in  the  program.  The  effects  are  assumed 
to  be  inversely  proportional  to  the  vertical  distance  from  the  wing, 
corresponding  to  thin  wing  theory. 

V ) ' ” ; ‘ . ■ i Tj  * ' * . 

Compressibility  effects  arise  only  in  the  calculation  of  velocities 
induced  by  the  line  vortices  used  to  model  the  circulation  effects.  To 
correct  for  compressibility,  all  dimensions  in  the  x direction  are  first 
multiplied  by  the  factor  1/  \Jl-Mz  where  M is  the  Mach  number  before  the 
calculations  for  induced  velocity  are  made.  (The  induced  velocity,  which 
is  circumferential  about  the  line  vortex,  would  be  Inversely  proportional 
to  the  distance  from  the  line  vortex  for  an  infinitely  long  vortex;  for  a 
finite  length  or  segment,  of  a line  vortex,  though,  it  Is  somewhat  less  and 
depends  on  the  distance  and  direction  from  the  point  to  the  ends  of  the 
vortex.)  The  vertical  and  sideways  induced  velocities  then  calculated  are 
correct  for  the  compressible  flow. 

. . i ■ • : ■ 

The  net  local  angle  of  attack  at  a given  point  1b  obtained  as  the 
sum  of  the  actual  angle  of  attack  (relative  to  an  airplane  reference  line) 
and  the  vertical  perturbation  velocities  obtained  from  the  three  sources 
mentioned  above  (since  U is  set  equal  to  1.0  and  the  perturbation  velocity 
in  the  x-direction  is  neglected  relative  to  U). 

>t.< v viljrtiji',.?0*T  afitzunq  ..bold-". -rib  e Yj 

The  net  local  angle  of  sideslip  is  obtained  as  the  sum  of  the  perturbation 
velocities  from  the  first  two  sources  mentioned  above  (there  not  being 
any  sidewash  effect  due  to  wing  thickness). 

Corrections  to  the  calculated  angles  of  attack  and  sideslip  are  made 
in  the  near  vicinity  of  the  wing  and  body  to  make  the  flow  approach 
tangency  to  the  wing  or  body  surface  as  the  distance  to  these  surfaces 
diminishes. 

The  inputs  to  the  program  and  corresponding  formats  are  given  below, 
with  a description  of  the  input  parameters  following: 

FORMAT 


1st  card: 

NX,  NY,  NZ,  NR 

415 

2nd  card: 

S,  DM,  B,  C4SWP,  CR,  XCR4,  TAU,  ZF 

8F10.4 

3rd  card: 

X(l),  1 - 1,  NX 

8F10.4 

' 7'h  Yj 

8F10.4 

’ n . ' r 'i.  7 * o n 

4th  card: 

Y(J),  J - 1,  NY 

• 

5th  card: 

Z(K),  K - 1,  NZ 

8F10.4 

6th  card: 

AM,  CL,  AL,  E,  DELRM1,  DELRM2 

6F10.4 

bi! 


3 


* 


* I 


(This  last  card  to  be  repeated  NR  times.) 


J 
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NX,  NY,  NZ  - number  of  x-values,  y- values,  and  z-values  respectively  In  a 
matrix  of  points  (including  all  possible  combinations)  at  each  of  which 
it  is  desired  to  calculate  the  local  angles  of  attack  and  sideslip.  NX, 

NY,  and  NZ  must  be  in  the  range  of  1 to  8. 

NR  - number  of  runs  (different  flight  conditions). 

S - theoretical  wing  area. 

DM  - body  diameter  at  the  location  of  the  wing  (should  be  equal  to  the 
horizontal  distance  between  the  exposed  wing  roots  as  is  usfed  for  wing- 
body  interference  effect  calculations). 

B - theoretical  wing  span. 

C4SWP  - quarter  chord  sweep  angle  of  the  theoretical  wing  in  degrees. 

CR  - theoretical  wing  root  chord. 

XCR4  - x-location  of  the  quarter-chord  point  of  the  theoretical  wing 
root  chord. 

TAU  - thickness-to-chord  ratio  of  the  wing. 

ZF  - z-locatlon  of  the  fuselage  (or  body)  centerline  measured  relative 
to  the  plane  of  the  wing  (positive,  above  the  wing). 

X(l),  Y(J),  Z(K)  - x,  y,  and  z-coordinates  of  the  matrix  of  points  at 
which  local  angles  of  attack  and  sideslip  are  to  be  calculated.  X-values 
are  referenced  according  to  the  value  of  XCR4  (more  positive,  farther  aft), 
Y-values  are  referenced  to  the  airplane  vertical  plane  of  symmetry 
(positive,  out  right  wing),  and  Z-values  are  referenced  to  the  plane  of 
the  wing  (positive,  above  the  wing).  X(I),  Y(J),  and  Z(K)  should  not  be 
such  that  there  are  points  in  the  plane  of  the  wing  (Z(K)  * 0.)  or  on 
the  fuselage  centerline  because  of  singularities  in  the  mathematical 
model  there.  Also,  Y(J)  must  be  positive.  (For  desired  points  on  the  left 
side,  input  the  corresponding  point  on  the  right  side.  The  results,  in 
terms  of  inboard  or  outboard  sidewash,  will  be  the  same.) 

AM  - Mach  number  at  the  flight  condition  in  question. 

CL  - lift  coefficient,  based  on  the  theoretical  wing  area  S,  at  the  flight 
condition  in  question. 

AL  - airplane  angle  of  attack  (relative  to  an  airplane  reference  line)  in 
degrees. 

K - n factor  described  above  equal  to  the  maximum  A (c^c)/cflV  (occurring 

at  the  three-quarter  semispan  point)  from  an  elliptical  (neglecting  wing-body 
Interference  effects)  lift  distribution  at  the  flight  condition  in  question 
(positive  if  lift  is  shifted  from. outboard  to  inboard). 

NOTE:  Effect  of  E is  small  so  only  need  to  input  if  there  is  a significant 
change  from  elliptical  loading  (not  counting  wing-body  interference) 


l 


7 
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I 
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DELRM1,  DELRM2  - increments  to  be  added  to  the  nominal  relative  magnitudes  of 
the  line  vortex  strengths  of  0.503  for  the  line  vortex  positioned  at 
0.135c  and  of  0.345  for  the  line  vortex  positioned  at  0. 25c, respectively, 
as  described  above.  (The  sum  of  these  two  relative  magnitudes  plus  that 
of  the  third  line  vortex  located  at  0.63c  (set  automatically)  equals  1.0.) 
This  defines  the  chordwise  distribution  of  lift. 

NOTE:  Since  effect  of  chordwise  lift  distribution  is  generally  small, 
only  need  to  input  for  DELRM1  and  DELRM2  if  there  is  q’jite  a large  change 
from  the  nominally  assumed  flat  plate  loading. 

For  all  inputs  of  lengths  or  areas  above,  the  same  unit  (or  this 
unit  squared)  should  be  used. 

Besides  the  input  on  data  cards,  certain  variables  must  be  defined 
with  statements  within  the  program  (before  they  are  used,  obviously). 

RF(I),  which  corresponds  to  the  body  radius  R(x)  at  the  x-location 
corresponding  to  X(I),  and  RPF(I),  which  corresponds  to  R/ (x) , must  be 
defined  for  all  NX  values  of  I.  (Most  often,  with  the  x-locations  in 
question  near  the  wing,  RF(I)  will  be  equal  to  0.5  x DM;  another  note  is 
that  in  the  case  of  fuselages  with  oblong  cross-sections,  the  horizontal 
dimension  is  usually  what  should  be  used  (at  least  with  regard  to  the  flow 
field  away  from  the  fuselage).) 

The  effective  fuselage  angle  of  attack  ALF  and  the  wing  angle  of  attack 
ALW,  which  can  vary  along  the  span,  must  also  be  defined.  ALF  is  typically 
defined  as  the  sum  of  AL  plus  any  incidence  of  the  fuselage  relative  to 
the  airplane  reference  line  (in  degrees)  minus  an  induced  angle  of  attack 
ALI  due  to  wing  downwash.  ALI  is  nominally  set  by  the  program  based  on 
the  expected  downwash  at  the  wing  for  the  given  lift  coefficient,  area, 
and  span.  The  wing  angle  of  attack  ALW,  as  presently  set  up  in  the  program, 
is  calculated  as  the  sum  of  AL  plus  EO,  the  wing  incidence  relative  to  the 
airplane  reference  line,  plus  a wing  twist  factor  ET  (positive  for  washing) 
times  the  relative  distance  out  the  semispan.  (EO  and  ET  are  in  degrees 
and  must  be  defined  in  the  program.) 

NOTE:  ALW  is  used  only  to  correct  the  flow  field  near  the  wing  and  does  not 
affect  other  parts  of  the  field;  thus  it  is  not  critically  important. 


FINAL  NOTES 

It  must  be  remembered,  that  even  though  it  is  done  on  a computer,  the 
flow  field  calculated  using  this  program  is  only  a rough  approximation. 
Factors  affecting  this  are  the  fact  that  the  geometry  is  simplified  from  the 
actual  case,  (usually,  slender  body  and  thin  wing  theory  are  used),  the  line 

vortices  which  are  used  to  generate  circulation  effects  are  finite  in  number 
and  the  lift  distribution  (including  wing-body  interference  effects)  is 
only  estimated;  also,  small  angle  approximations  are  used  and  the  change 
in  longitudinal  velocity  relative  to  the  free-stream  velocity  is  neglected 
in  the  calculations  of  downwash  and  sidewash  angles.  In  addition, 
influences  from  other  parts  of  the  airplane,  such  as  tail  surfaces,  are 
neglected.  A typical  error  might  be  about  0.5  degrees  (although  in  certain 
places  it  could  be  more),  as  indicated  by  comparison  of  the  results  in 
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Appendix  B,  which  are  for  the  A-7  airplane,  with  wind  tunnel  data  reported 
in  reference  (c).  The  wind  tunnel  results  themselves,  at  least  as  presented 
in  that  reference,  are  subject  to  errors  of  about  the  same  magnitude. 

t 

It  is  possible,  using  the  results  from  this  program  as  a first  order 
correction  to  free  stream,  to  further  refine  the  estimate  of  the  flow  field. 
As  an  example,  the  effect  of  a pod  can  be  estimated  by  placing  it  in  the 
flow  field  as  determined  by  this  program  and  then  using  slender  body  theory. 
Also,  since  the  separate  effects  of  the  body,  wing-body  circulation,  and 
wing  thickness  on  the  flow  disturbance  are  indicated  in  the  output,  local 
corrections  for  any  of  these  can  be  made. 
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FIGURE  1.  ARRANGEMENT  OF  LINE  VORTICES  ON  WING 
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,r  - Fr0U-(y/s)' 


_______ 

rr  - <>|l-<r/s)2~a)] 


Jlliptical  distribution 

r - r0\ji-(y/s)2 


r * body  radius 


S ■ wing  tip 


Heavy  line  indicates  modified  lift  distribution  to  account  for 
wing-body  Interference. 


FIGURE  2.  SPANWISE  LIFT  DISTRIBUTION  CHANGE  (AS  MODELED) 
DUE  TO  WING-BODY  INTERFERENCE 
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APPENDIX  A 
LISTING  OF  PROGRAM  WBFLFD 


A-l 


♦CO 


t *1 . ✓ ( 1 .-*S*r. d t J*b</HT  1 A ,./KS*HS/* . ) ♦ c 1 . *H)  > 
IFtuS.HF.O.)  **l.!»71/RS»H/Cl.*H)-St>ftTO.-BS*RS) 

00  »0  t«1»nx 

wri*.  *»*om 

I t 1 ■ 0 • 

1 0 ■ - ] , 

£y__-  . 

oo  so  Miiiin 

YT (M)«v#.0S*P 

VtfPiJ  TvT7ST".Sg*Vfi • — • T”T~ 

*'»<•*»  I ) s*C»4-.  1 tt»*CK*Y<i  I*)  * <TLC4*.?3/H#CC> 

VH  I s )c>f.  ) • TIC* 

' AtTJu  « i)  «UU<i».to/H*CC) 

uo  i. lo 

l»0  Sb  »‘*li3 

If  (vi.( [l-L i .iT.b^/'Zr.J  x:»ni-L*N) »xft<  l«i*t  »N) 

I/O  rM>  f*n*i»mi/ 

Hr. Alt  m»AF»Ci  ,»tit  *0fct»Ml  ,OFLLMi* 
mm  I 1 

“«<■»>«.  I sy-pFi.v*  I -Dti.s*a 


►MINT  *0 


• '*N  • '• 


so  1-uWM.T  <////. SX»«mA Ch  FO.  LIFT  COfeM.  ALPHA  tOFi*)  . f.  •♦SOM 
Of LTA  T | CL *C) LOCAL ) /C  AV«4Xt»Ww(n  RNTP)  P»M3) 

J 'JiTt.  OF*?s|  H.LlPTlLAL  JNHU  SHIFTm./'H* **CHH'Ust  l! 

♦Ff  HJST.«/)UX«*AT  tl3NC.  #i»H» 

y«lUX  •>*•«  A"»M.  *AL  *t  •H*' ( 1 1 tb"f<  ) *•*"  t j)  . 

~s  -rr:vTfTTS  rrr~;7  |,1*  jyx»iH.3> 

n|At‘iu<T(|,«/.M»»r)  . . 

CL»U.»«ft«,»l/w#|>M/S»AL  . ......  *■  v.,. 

*1  I *rf>  ♦ •< 

.»»  F **|  -At,  I . 

• It)  H>  '<*1  ♦ | o - l L 1 

rs«Yo*”)  /•♦•/»# 

IF  r.<*(i»)>F«liO»<>>kT(U»Y<tTSl  v ■ » 

"Vr  i «hsi  <.-<  <m |/u  <sv»  t 1 1 .-’/sm'SI-a) ttys/H*. ) 

IF  (I  •/./' (0« • A/  0» YS »f>T« »S)  r.H0J)aaH<H)*.s*r*C*MMF.*c3»YS> 

A-* 


APPENDIX  A 

PHOUPAH  PMFLFul  INPUT  tOUtPtil) 

DIMENSION  XCH)  *Y|b) »II8)  ,PM(S>  ,RF(M>  .PPFI8>  fYT.ClO)  ♦ VMlO)  * 

•»£Ap.  h.MfMV.WtNW  - , h 

FUPmaI  (Alb) 

HEAD  1S»S *L)M(p tCASwPfCR* xrRA *TAU*ZF  * .>  ; 

*£!«*»  I * ■■  - 

PHINT  <»0 

FOPmaT ()Hl»SXf*MlNG  APtA#,3Xt*H0DY  Ol AM**3X**WlN0  SPAN#f3X 
VC/H  St.t>p,(OFB)  FOOT  CHORD  X AT  CP/»  THKSFSS  J» 

3 **••!  nr  HOOV  CTRL INF */) 

H«InT  /*s»«$«UM»b«C‘»S*r>*CR»xCR*  tTAUtZF 

_ti^ii*Tjj>F_ip..?,F  !♦.  lyFUy  ItS  2 

mr  •»!*  1“.  <*!!)«  ) 

NfcAU  1S.(Y(J) ,J«1»NT) 

Ht All  |S, (7 <K) ,K«1 »N7)  , 

C«s/H  ..  . ...  , V»W  J'-iy, 

CCs*. *<0-0 

MS«0*/P  . 

lLC*=  f af(C‘»<>wm/;>7.2ss) 

ILLF.«TLCA*.b/H*CC 

TUIF.«T|.rA-I.S/M*CC 

»**An  I « ' ' . " 

IF  J ritF  .liT  . 0 , ) 0»CCM-PS*fC)/iCP*.S*<F'-nM)*TLC<t-<,7S*kS/*.) 


xotift*: 


BESfAVAl 


Fi>**ATl///,7X.»Y«tt)X.«X*.MX.*7»»SX»MU»M#  LOCAL  ’lOfG) 
^ HiH’l 


BETA  LOCAL 


ALOHA  uue 


<POS.»  QUTHOAHU7 10U 


Art  CIHC  TO  HOny  TO  BOOT  ' Wiw»  TH*fNfSS»/» 
i)0  I'tu)  .l«l«MV  


YYS*Y ( J) *VT |M) 
DO  1 fill  ►-«!  ♦ 3 


X«1*<x<I)-X0<m«1iM l/BTA 

xxl«xx 


uUSI^T  (XXMXX)»YY»YY«?  (m*Z(K)) 


FF*i,h<m)»H m(X  ) •^.bSV 


W*r*K>FC*<  IYI.|)-YT<M-1J  )/SOPTIXX*XXMY|J)-YT|M-l>  >• 

>-(Ylw)  *VT  | ) /S'JHT  IX>*XX»«<Y  IQ)  *YT  lv»i)  ) 


H<?«YY*YV*7I*J*7IM 

h><'«yy«;*vvs*7ik)*7ik» 


IF IM.LF .V)  Vw«Vw*FF*2  <K » • | |XX/0»-XXi/0Sl) /MS2- t XX/O-XXl/Ol I /**) 

IF  |fi#>.»  ,f  I r.c,  Tn  100 


'■S<*<XXl/l>M*lTn 


IK)  «|  (»Xl/USl«l  •>/**?- (XXl/OWl  ,)/Pk!) 
COKTlMiF  • 


Bit  aT»Pf  n J A/JP'VT^vfrrULFAkf  ( 1 ) • (Y  iji  Vy  ( J)  -S.rF*t  7F>  /Li  ><M-F  | (TTk? 
VF«|hmf  (1)*Y(J)*ST.C9o-AlF«HF  ll)*?.*7ZF*Y(J)/PiO«HF  I |l 


RF.P//..OH.YIO) «|  T«UHS*«)  00  TO  ?Qtv 


IF  IT  |J)  . 

XC«IX  m-XCo-.C.w/^.-y  (J|«TLLt  l/CH 

I f* I XI  «|  » .U.»nw.<C«l-T«l  »)  »0  TO  •■«  

r^u7.-s(>Tnr;r  t . • *•.•<>  / it  a i » •••l***»*a»«  «/  ifj  i/r**i*T**in«At» » 

#»1h« I am*TAM/7 In  I #LH* (II  (l7.MdAXC-Ae*bA>*xC««Bk0M*xr.-SO.I  j)  *xc* 

) IF  IM  .LF.O.)  00  TO  *ll> 

MmF ■M^'TfHiO/HF  III  

rrU't'T  .TTTTTVTTNiri  YiJl-Yfjl  /‘o«;.Lt  «.iik  Hk«-ALl**.*i*oF-i ,» •!**»•  *l 
IFIMwF.i  T,1  .‘-.X'lO.Vl  J)*Ym)/P7.'1T..!>|  VY*F.  0*vw»  (DMF-I  • ) 
i m.iy.  »AL 
/ «v • ♦ vT  " 
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APPENDIX  B 
SAMPLE  OUTPUT  FROM  PROGRAM  WBFLFD 


KLTAIICU*C»LOCM.I/C  AV  (Mil)  »*tJ>  (M*$X«CL.  **•.  of 

flLMtCM.  !•*««.•  !*«tO  SHl»TI  CMMXSF  LlF?  OI*T. 

.m  .1^ 


•»J 


